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Abstract 

(N 

^ ' We give the basic definitions and some theoretical results about hyperdeter- 

. minants, introduced by A. Cayley in 1845. We prove integrability (understood 

as 4d-consistency) of a nonlinear difference equation defined by the 2x2x2- hy- 
' perdeterminant. This result gives rise to the following hypothesis: the difference 

ff^ ■ equations defined by hyperdeterminants of any size are integrable. 

I We show that this hypothesis already fails in the case of the 2x2x2x2- 

hyper determinant. 

> 

X ■ 1 Introduction 



Discrete integrable equations have become a very vivid topic in the last decade. A 
number of important results on the classification of different classes of such equations, 
based on the notion of consistency [3] , were obtained in [H |l7j (cf . also references to 
earlier publications given there). As a rule, discrete equations describe relations on the 
scalar field variables G C associated with the points of a lattice with vertices 

at integer points in the n-dimensional space = {(xi, . . . ,Xn)\xs G M}. If we take 
the elementary cubic cell Kn = {{ii, . . . ,in) \ is ^ {0, 1}} of this lattice and the field 
variables fii...i„ associated to its 2" vertices, an n-dimensional discrete system of the 
type considered here is given by an equation of the form 

Qn(f) = 0. (1) 
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Hereafter we use the short notation f for the set (/00...O5 ■ • ■ , /11...1) of all these 2" 
variables. For the other elementary cubic cells of the equation is the same, after 
shifting the indices of f suitably. 

The equations mostly investigated so far [H [21 [T7j were supposed to have the fol- 
lowing properties: 

1) Quasilinearity. Equation ([1]) is afiine linear w.r.t. every fi^i2...ini i-G- Q has 
degree 1 in any of its four variables. 

2) Symmetry. Equation ([T]) should be invariant w.r.t. the symmetry group of 
elementary cubic cell Kn or its suitably chosen subgroup. 

On the other hand a number of interesting discrete equations which do not enjoy 
one or both of these properties has been found. In this publication we investigate an 
important class of symmetric discrete equations which do not have the quasilinearity 
property and are given by the equations Hn{foo...o, ■ ■ ■ , = 0, were Hn denotes the 

n-dimensional hyperdeterminant of the corresponding n-index array (/00...O! • • • ? 
We give the precise definition of hyperdeterminants in Section [2l In the simplest two- 
dimensional case of the 2x2 matrix the hyperdeterminant is just the usual determinant: 

H2{i) = /oo/ii - /oi/io. (2) 
The next nontrivial case is the 3-dimensional 2x2x2- hyperdeterminant: 

-f^3(f) = /ni/000 + /loo/oii + /loi/oio + /iio/ooi 
~2/iii/iio/ooi/ooo ~ 2/111/101/010/000 
~2/iii/ioo/oii/ooo ~ 2/110/101/010/001 (3) 
~2/iio/ioo/oii/ooi ~ 2/101/100/011/010 
+4/111/100/010/001 + 4/110/101/011/000- 



/oil /ill 




Figure 1: Square K2. Figure 2: Cube K^. 

The corresponding elementary cells K2, K3 and the field variables associated with 
the vertices are shown on Figures [H El 

The general definition of hyperdeterminants was given by A. Cayley [7j, who also 
gave the explicit form ([3]) of the first nontrivial 2x2x2- hyperdeterminant. In the last 
decades, following the modern and much more general approach of ^-discriminants [9] , 
the theory of hyperdeterminants found important applications in quantum informatics 

biomathematics [S], numerical analysis and data analysis [6] as well as other fields. 

As one can easily see, the expressions (3.7) in [TB] and (6.11) in [12], describing 
some discrete integrable equations, are nothing but the classical Cayley's 2x2x2- 
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hyperdeterminant ([3]) . We prove below in Section [3] that ([3]) is also integrable in the 
sense of (n + 1) -dimensional consistency [3]: 

An n- dimensional discrete equation (J\) is called consistent, if it may be imposed in 
a consistent way on all n-dimensional faces of a (n + 1)- dimensional cube. 

We give the accurate formulation of this general consistency principle for the case 
of non-quasilinear expressions similar to (j3]) in Section [31 For the two-dimensional de- 
terminant ([2]) (which is quasilinear) consistency can be established by a trivial compu- 
tation; the equation H2{foo, fii, foi, /lo) = is obviously linearized by the substitution 
fij = exp fij . Using a result on Principal Minor Assignment Problem proved in \1AI 
we establish 4(i-consistency of the 2x2x2- hyperdeterminant ([3]) in Section [3], cf. 
Theorem [2] below for the precise formulation. 

This result gives rise to the following Conjecture: the difference equations defined 
by hyperdeterminants of any size are integrable in the sense of (?t, + l)-dimensional con- 
sistency. Nevertheless as we show in Section HI this Conjecture fails already in the case 
of the 2x2x2x2 - hyperdeterminant. The computation of this Ad - hyperdeterminant 
turns out to be highly nontrivial (compared to the relatively simple expressions ([2]), 
(IHl)) and was completed only recently [13]. We report in Section H] a more straightfor- 
ward and simpler computation of the same hyperdeterminant with the free symbolic 
computation program FoRM [18j. The size of this hyperdeterminant (2 894 276 terms, 
total degree 24, degree 9 w.r.t. each of the field variables) implies that checking its 5d- 
consistency can be done only numerically, using high precision computation of roots 
of respective polynomial equations on the 4c?-faces of the 5-dimensional cube K^,. This 
was done using two different computer algebra systems Reduce [T^ and Singular 
[20j . As our computations have shown (cf. their description in Section Hj), the Ad - 
hyperdeterminantal equation H^li) = is not 5(i-consistent. This non-integrability 
result should be investigated further since recent examples [11] show that consistency 
is not the only possible definition for discrete integrability. 

2 The definition of hyperdeterminants and its vari- 
ations 

The remarkable definition of hyperdeterminants given by A. Cayley in 1845 [7] and still 
used today |9j describes the condition of singularity of an appropriate hypersurface. 
Let A = (aj^i2...j^) be a hypermatrix (an array with r indices) with = 0, . . . , n^. The 
polylinear form 

7-r _ (1) (r) 

h---ir 

defines a hypersurface f/ = in CP"^ x . . . x CP"''. Here xp"* denote the homogeneous 
coordinates in the respective complex projective space CP"^*. This hypersurface is 
singular, i.e. has at least one point where the condition of smoothness is not satisfied 
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iff tlie following set of {rii + 1) ■ . . . ■ {rir + 1) equations 
<Vs = l,...,r, \fk=l,...,ns, 




dU 







(4) 



has a nontrivial solution E CP"''. As one can show (cf. [9]), if a certain condition ([5]) 
on the dimensions of the array A is satisfied, elimination of the variables xf^^ from (jl]) 
results in a single polynomial equation in the array elements aj^ij - ir- Hr{A) = 0. This 
polynomial is irreducible and enjoys practically the same symmetry properties as the 
usual determinant of a square matrix. Following Cayley this polynomial Hr{A) (defined 
uniquely up to a constant factor) is called the hyperdeterminant of the array A. The 
necessary and sufficient condition of existence of a single polynomial condition Hr{A) = 
for the hypersurface ?7 = to be singular, i.e. the condition for the corresponding 
hyperdeterminant of A to be correctly defined, is as follows: 



In particular, if r = 2, so for usual (ni+l) x (ri2+l)-matrices, this condition implies ni = 
n2, and in this case the hyperdeterminant H2 coincides with the classical determinant 
of the matrix Ai-^i^. Note that for a given set {rii, . . . , rir} of array dimensions one says 
that we have the corresponding (ni + 1) x . . . x (^,.+1) - hyperdeterminant since the array 
indices range from to Uk- The hyperdeterminant is SL{C, rii + l) x ■ ■ ■ x SL{C, rir + l)- 
invariant, which means that if one adds to one slice A^^p = {(aiiia - ir) with fixed ik = 
p} another parallel slice A^^g, q p, multiplied by some constant A, the value of if^ is 
unchanged; swapping the slices A^^p, A^^q either leave again invariant or changes its 
sign depending on the parity of the dimensions rij; finally, multiplication of a slice A^ p 
with a constant A results in multiplication of the hyperdeterminant by an appropriate 
power of A. Hr is also invariant w.r.t. the transposition of any two indices i/, im of the 
hypermatrix A = (aj^jj-v)- 

As we have stated in the introduction, the first nontrivial 2x2x2- hyperde- 
terminant ([2]) was computed by A. Cayley himself [7|. Amazingly enough, already 
the next step, computation of the 2x2x2x2 - hyperdeterminant is very difficult. 
The problem of computation of an explicit polynomial expression for this case was 
proposed by I. M. Gel'fand in his Fall 2005 research seminar at Rutgers University. 
The monomial expansion of the 2x2x2x2 - hyperdeterminant is related to some 
combinatorial problems, and was done (using an inductive algorithm of L. Schlafii [_8j) 
for the first time in [TS] , using a dedicated C code; this computation required a serious 
programming effort since the standard computer algebra systems like Maple can not 
cope with the intermediate large expressions. The resulting polynomial expression for 
the 2x2x2x2 - hyperdeterminant has 2 894 276 terms, total degree 24, and has 
degree 9 w.r.t. each of the array entries aj^jjisM- The size of this expression in usual 
text format is around 200 megabytes. 

In October 2007 we re-checked this computation of the 2x2x2x2 - hyperdeter- 
minant using a free symbolic computation program Form [18] and the same inductive 




(5) 
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algorithm of L. Schlafli [8]. The computation required 8 hours on a 3GHz processor 
of one of the SHARCNET nodes (www.sharcnet.ca) and some 800 Mb of temporary 
disk storage. Due to the efficient design of Form one had no need to write any special 
low-level code, the standard Form routines are completely sufficient. The obtained 
expression can be saved either in text format (around 200 Mb) or in the internal binary 
Form format; both can be used as input to other FORM runs and require only around 
1 minute to be read into a session. Moreover, Form could cope with a straightforward 
check of the invariance of the obtained hyperdeterminant w.r.t. addition of one slice 
Afc p multiplied by a symbolic constant A to another parallel slice Ak,q, q ^ P- 

This more challenging computation still used the usual Form routines and required 
some 10 hours of CPU time and around 200 Gb of temporary disk storage (the size 
of the intermediate expression reached ~ 800 000 000 terms!). This may suggest that 
some other hyperdeterminants of higher size may also be computed using available 
software and hardware. 

Amazingly enough, if one reads the original Cayley papers [7j, then a relatively 
small expression with some 340 terms which Cayley describes as the 2x2x2x2- 
hyperdeterminant can be found! A straightforward check shows that this expression 
enjoys the invariance properties for hyperdeterminants stated above. On the other 
hand. Proposition 1.6 in ^ p. 447] states that if a polynomial in the entries aj^jj-' ir of 
a hypermatrix A has these invariance properties and meets some extra weak condition 
on the stars of the monomial powers then it should be divisible by the respective 
hyperdeterminant of A. The Cayley's expression does satisfy the necessary invariance 
conditions; only a few terms do not meet the required extra condition on the stars of 
the monomial powers, so this statement of Cayley on the explicit form of the 2x2x2x2 
- hyperdeterminant is wrong. This was already remarked by L. Schlafli [8] who gave 
an inductive algorithm for the computation of hyperdeterminants. Unfortunately, as 
shown in [10] , the Schlafli algorithm works only for very special hyperdeterminants, in 
particular the only hyperdeterminants with Ui = . . . = = 1 which it can compute 
are precisely the 2x2x2- and 2x2x2x2 - hyperdeterminants. Already for the 
2x2x2x2x2 - hyperdeterminant there seems to be no better way to compute 
the explicit expression other than the elimination procedure given in the definition of 
hyperdeterminants. 

Based on the above symmetry discussion one could adopt other definitions for 
hyperdeterminants. Historically this resulted in a few other definitions of hyperdeter- 
minants as invariant expressions. Many of them have a much simpler form than the 
definition we apply. A review of various definitions can be found in |12j . 
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3 The Cayley 2x2x2- hyperdeterminant as a dis- 
crete integrable system and the Principal Minor 
Assignment Problem 

The integrability definition applied in this pubhcation is based on the requirement 
of consistency which is described in this section in detail for the case of the 3d - 
hyperdeterminant ([3]). Suppose we have a 4d-cuhe shown on figure H] with field 
values fijku hj,k,l G {0,1}. One should impose the formula ([3]) on every 3(i-face 
of K4, by fixing one of the indices i,j,k,l, and making it for the faces which we 
will call below "initial faces", or respectively 1 for the faces which we will call "final 
faces". Further one needs to fix some mapping from the initial "standard" Sd-cube 
shown on figure [3] (with the vertices labelled fijk) onto every one of the eight 3d- 
faces (for example {fnki} on {x2 = 1}). Due to the symmetry properties of the 3d 
- hyperdeterminant ([3]) this can be done, for example, using the trivial lexicographic 
correspondence of the type fijk ^ fujk- The initial data are some arbitrary complex 
values of the field variables fijki assigned to the vertices shown on figure H] by black 
circles. Then, using the equation imposed on the initial faces we can find the values 
of the field variables for the last (8**^) vertex of the respective initial face, such vertices 
are shown on figures [3l H] as white circles. Obviously since the equations ([3]) for these 
values are quadratic we obtain two possible values. On the next step we impose ([3]) 
to hold on the 4 final faces using the initial data and the values found on the previous 
step. From each of the 4 final 3(i-faces we again obtain possible values for the final 
vertex /mi shown on figure H] by a small white box. For generic initial data we have 
on each final face 2^ = 8 different choices for the intermediate values (white circles) so 
we find in principle 8 ■ 2 = 16 different possible values for /mi from each final face. 
How many of them coincide among the 4 final faces? 



As we prove below, 8 of every 16 values for /mi found for each of the final 3(i-faces 
are common. All other 8 ■ 4 = 32 are in general not shared between the final faces as 
it has been confirmed by numerical examples. 

This result should be considered as the proof of consistency for face formula ([3]); 
our considerations below are based on a remarkable result proved in |14j . 




Figure 3: Cube 



Figure 4: Cube K4. 
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Let us first formulate tlie necessary definitions and results of [H]. Suppose we have 
a real symmetric nxn-matrix M = [rriij). Its principal minors form a vector of length 
2"^ with entries indexed by subsets / of the set {1,2, ... ,n}. Namely, Mj denotes the 
minor of M whose rows and columns are indexed by I. This includes the OxO-minor 
M0 = 1. The famous Principal Minor Assignment Problem considers the description of 
a suitable complete set of algebraic relations among the minors of a generic symmetric 
nxn-matrix M = [rriij). The first observation (formula (2)) of [H] consists in the fact 
that for principal minors of a 3 x 3 symmetric matrix one has the following relation: 

^0^123 + ^1 ^23 + ^1^13 + MlMf2 + 4M0M12M13M23 + 4M1M2M3M123 

-2M0M1M23M123 - 2M0M2M13M123 - 2Mti,M3MuMu3 - 2M1M2M13M23 
-2M1M3M12M23 - 2M2M3M12M13 = 0. 

This obviously gives us the Cayley's hyperdeterminant ([3]) if we interpret every minor 
Mj as the field variable /iiiaig with = 1 if s G / and = otherwise, for example 
= /loi- For the initial vertex we have /ooo = ^0 = 1- For symmetric matrices M 
of larger size the 2x2x2- hyperdeterminantal relations are also fulfilled for "shifted" 
principal minors, in our terminology this means that the hyperdeterminant ([3]) vanishes 
on every 3(i-face of the n-dimensional hypercube with field variables equal to the 

principal minors Mj such that is = 1 ii s & I and is = otherwise. As a remarkable 
fact (not necessary to us) we mention that for d > 2 on any c?- dimensional face of 
this n-dimensional hypercube the corresponding d-dimensional hyperdeterminant also 
vanishes: 

Theorem 1 ([14]) Let M = {rriij) be a symmetric nxn matrix. Then the vector M^, 
of all principal minors of M is a common zero of all the hyperdeterminants of formats 
from 2x2x2 up to 2x2x — - x2 . 

n terms 

For n > 2, the entries of the symmetric matrix M = {rriij) are determined up to sign 
by their principal minors of size 1x1 and 2x2 since ma = Mi, mfj = MiMj — MijM^. 
As we see, the value /m = M123 for the final vertex is not defined uniquely from 
the relation ([3]); the choice of one of the two possible values corresponds to a suitable 
choice of the signs for the 3 off-diagonal elements mi2, rrii^, m23. These 2^ = 8 sign 
combinations give only two different values for M123, because changing simultaneously 
the signs of a certain row and the column symmetric to it we do not change the principal 
minors of the matrix, so we can fix the signs of mi2 and mi3. 

If we take a 4 x 4 symmetric matrix, all its elements will be defined up to the 
sign of the off-diagonal elements by the principal minors of size 1x1 and 2x2; this 
corresponds precisely to the choice of the initial data (black circles) on figure HI So 
provided we choose the initial data fijki with /oooo = 1, we can find a corresponding 
4x4 symmetric matrix with fixed 1x1 and 2x2 minors; the off-diagonal elements of 
the matrix are fixed up to the signs. So, fixing 3 signs of the 6 off-diagonal elements 
(changing simultaneously the signs of a certain row and the column symmetric to it) 
we have only 3 essential sign choices; as the numerical examples show, they really give 
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2^ = 8 different minors M1234 = det M. Precisely these 8 values should give the 8 
values for /mi that should coincide after computation of 16 possible candidates for 
/nil from each of the 4 final faces from the initial data. 

If we choose initial data with /qooo 7^ 1; /oooo 7^ 0, homogeneity of the face equations 
([3]) allows us to reduce this situation to the case /oooo = 1 considered above. Thus the 
following statement has been proved: 

Theorem 2 Let some generic initial data {/oooo, /looo, /oioo, /ooio, /oooi, /iioo, /loio, 
/looi; /oiiO) /oloi; /ooli } on t/ic cn6e 6e given. After computation of the two possi- 
ble values for each of the intermediate vertices {/mo, /iioij /loii? /oiii} from the face 
relations ^ on the respective 3-dimensional initial faces, among the sets of 16 possible 
values of /mi for each if the 4 final faces, subsets of 8 values coincide for all of them. 
They are equal to the 8 possible values of det M for the symmetric matrices having 
given 1x1 and 2x2 principal minors corresponding to the initial data for fijki. 

4 The next step: the 2x2x2x2 - hyperdeterminant 
and its 5(i-inconsistency 

The definition of hypothetical 5(i-consistency for the 2x2x2x2 - hyperdeterminant 
which has degree 9 w.r.t. each of its 16 variables is easily formulated along the lines given 
in the previous Section. So in this case not only the size of the face equations, but also 
the number of possible choices of the intermediate values fijkim with i+j + k + l+m = 4 
for the computation of the final /nm from each of the 5 final faces of the 5-dimensional 
hypercube is dramatically increased. 

The strategy of numerical checking the hypothesis of 5(i-consistency adopted by us 
involved the following steps. 

1) We assign some random integer values for the initial data fijkim with i + j + k + 
I + m < 3 and use Form to substitute them into the expressions of the 2x2x2x2 
- hyperdeterminant on the 5 initial 4(i-faces obtaining univariate polynomials for each 
of the intermediate /oim, . . . , /imo and output the resulting expressions into a text 
file for further processing by REDUCE and independently by Singular. For initial 
data being random integers in the range [1, 100], the obtained equations have integer 
coefficients with approximately 40 decimal digits. 

2) The same substitution of the initial data into the final faces is performed with 
Form, resulting in much larger multivariate polynomials for the intermediate values 
/oim, . . . , /imo and the final /nm. These polynomials (each has the size of ~360 Kb) 
are output into a text file for further processing by Reduce and independently by 
Singular. 

3) The 5 univariate polynomial equations for the intermediate /oim, . . . , /imo are 
solved with a guaranteed precision of 150 digits. 

4) For each one of the 9^ combinations of the 9 complex roots the following com- 
putation is performed. 
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a) The set of complex roots for /oiiii, /iiiio is replaced in the 5 final face 
relations obtained on step 2, which makes them univariate polynomials for /nm with 
complex rounded coefficients. 

Starting with a guaranteed precision of 20 digits: 

b) One of the obtained polynomials is solved for /nm and 

c) successively the other 4 polynomials are solved for /nm as long as there is a 
non-empty approximate intersection (with a definite relative tolerance, see below) of 
the sets of roots for the /nm for all the polynomials solved so far. 

d) If all five face relations have at least one common approximate solution then 
execution continues with step b) with twice as many guaranteed precise digits, up to 
a maximum of 80 digits. This was never necessary. 

To increase safety, two complex values u, v were only considered NOT to be ap- 
proximately equal if for p precise digits the difference u — v differed significantly from 
zero, more precisely if \u — v\/\u\ > 10"^/^. 

The computation performed on the nodes of SHARCNET has shown that no 
equal values for fmu are obtained from the 5 final faces. The details of the com- 
putation and the code used can be obtained from the authors or downloaded from 
http : / /lie .math . brocku . ca/twolf /paper s/TsWo2008/. 

This results in the conclusion that: 

The 2x2x2x2 - hyperdeterminant is not bd- consistent. 

The safety of our numerical inconsistency result is increased by the fact that it is ob- 
tained by two completely different computer algebra systems. The Singular package 
is written in C using an arbitrary precision C library for their numerical computations 
whereas Reduce uses a long number arithmetic implemented in Lisp. 



We have also checked that the relatively small expression given by Cayley [7j for 
the 2x2x2x2 - hyperdeterminant is also 5(i-inconsistent. 
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